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Abstract 

Effects of the anisotropy of Cooper pairs in spin-triplet ferromagnetic su- 
perconductors are investigated on the basis of the Ginzburg-Landau theory. A 
special attention is paid to the triggering of the superconducting state by the 
ferromagnetic order. The ground states of these superconductors are outlined 
and discussed. The idea about a possible coexistence of ferromagnetism and 
spin-triplet superconductivity in neutron stars is introduced. 

1. Introduction 

The cold and dense matter in the interior of neutron stars undergoes phase transitions 
to both superfluid and superconducting states (see, e.g., Refs. [U 12 El- Both the 
superfluidity due to neutron Cooper pairing and the superconductivity due to proton 
Cooper pairing may be of unconventional spin-triplet (p-wave) type [h. The latter 
is known from the theory of 3 He liquids |H EJ E] as well as from the theory of heavy 
fermion [HE]) and high-temperature ( EH HU H2] ) superconductors. 

The possible superconducting phases in unconventional superconductors are described 
in the framework of the general Ginzburg-Landau (GL) effective free energy func- 
tional ^T] with the help of the symmetry group theory. Thus a variety of possible su- 
perconducting orderings were predicted for different crystal structures O EH UH EES] ■ 
A detailed thermodynamic analysis [5J E] of the homogeneous (Meissner) phases and 
a renormalization group investigation 9] of the superconducting phase transition up 
to the two-loop approximation fl] were also performed. 

Recent experiments ^Zj at low temperatures (T ~ 1 K) and high pressure (T ~ 1 
GPa) demonstrated the existence of spin triplet superconducting states in the metallic 
compound UGe 2 . This superconductivity is triggered by the spontaneous magnetiza- 
tion (M) of the ferromagnetic phase (M-trigger effect JEj)- The ferromagnetic order 
exists at much higher temperatures and coexists with the superconducting phase in the 
whole domain of existence of the latter below T ~ 1 K; see also experiments published 
in Refs. JIH] • Moreover, the same phenomenon of existence of superconductivity at low 
temperatures and high pressure in the domain of the (T, P) phase diagram where the 
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ferromagnetic order is present has been observed in ZrZn and URhGe, too; for details, 
see, e.g. Ref. (221 El EH EE] • Note, that the superconductivity in the metallic com- 
pounds mentioned above, always coexists with the ferromagnetic order and is enhanced 
by the latter. Besides, in these systems the superconductivity seems to arise from the 
same electrons that create the band magnetism. 

A similar phenomenon of coexistence of ferromagnetism and spin-triplet superconduc- 
tivity may exist also in neutron stars. In this case the Cooper pairing of fermions 
(protons) will be triggered by the spontaneous magnetic moment of the same proton 
subsystem of the nuclear star matter. The basic features of these phenomena can 
be described within an extension of the GL theory. The results can be applied to a 
new and interesting problem of coexistence of superconductivity and ferromagnetism 
in neutron stars. The coexistence of spin-triplet superconductivity and ferromagnetism 
may explain the large magnetic field of these stars. To our best knowledge the problem 
of a possible coexistence of ferromagnetism and superconductivity in neutron stars and 
asprophysics objects is introduced for the first time in our present report. 

The theory allows the description of various types of phase transitions and multicritical 
points [TTIESI- Following notations in Ref. [IE], here we summarize previous studies [221 
123} l2~H HHj and present new results on the effect of the Cooper pair anisotropy on the 
phase diagram of spin-triplet ferromagnetic superconductors [TH] . 

Our consideration is focussed on the ground state, namely, we are interested in uniform 
phases, where the order parameters (the superconducting order parameter ip and the 
magnetization vector M = {Mj,j = 1, 2, 3}, do not depend on the spatial vector x G V 
{V is the volume of the system). 

2. Ginzburg-Landau free energy 

Consider the GL free energy F(ip, M) = Vf(ip, M) , where the free energy density 
f(ip,M) (for short hereafter called "free energy") of a spin-triplet ferromagnetic su- 
perconductor is given by [THj 

f^M) = a s H 2 + ||^| 4 + y|^ 2 | 2 + |X:^| 4 + G / M2 + y M4 C 1 ) 

j'=i 

+i^M.(ip x $*) + SM 2 \ij\ 2 . 

In Eq. (1), tp = {ipj]j = 1, 2, 3} is the three-dimensional complex vector = ip'^+iip") 
describing the unconventional (spin-triplet) superconducting order and B = (H + 
4ttM) = V x A is the magnetic induction; H = {Hj;j = 1,2,3} is the external 
magnetic field, A = {A,; j = 1,2,3} is the magnetic vector potential (V. A = 0). In 
Eq. (1), b s > 0, bf > 0, a/ = aif(T — Tf) is given by the positive material parameter af 
and the ferromagnetic critical temperature Tf corresponding to a simple superconduc- 
tor (M = 0), and a s = a s (T —T s ), where a s is another positive material parameter and 
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T s is the critical temperature of a standard second order phase transition which may 
occur at \H\ = Ai = 0; Ai = \M\. The parameter u s describes the anisotropy of the 
spin-triplet Cooper pair whereas the crystal anisotropy is described by the parameter 

The two orders - the magnetization vector M = {Mj} and xjj = {A/}, interact through 
the last two terms in (1). The 70— term [21] ensures the triggering of the supercon- 
ductivity by the ferromagnetic order (70 > 0) whereas the 5— term makes the model 
more realistic in the strong coupling limit (2D]- Both ?/>M-interaction terms included 
in (1) are important for a correct description of the temperature-pressure (T, P) phase 
diagram of the ferromagnetic superconductor ^Hj- In general, the parameter 5 for 
ferromagnetic superconductors may take both positive and negative values. 

The values of the material parameters (T s ,Tf,a s ,af,b s ,u s ,v s ,bf,Kj,j and 5) depend 
on the choice of the concrete substance and on intensive thermodynamic parameters, 
such as the temperature T and the pressure P. One may assume that the general 
form (1) of the free energy may describe the general features of the uniform orders in 
neutron stars provided one makes a suitable choice of parameters (T s , Tf, a s ,...). 

As we are interested in the ground state properties, we set the external magnetic field 
equal to zero (H = 0). Besides, we emphasize that the magnetization vector M may 
produce vortex superconducting phase [TJ [TT] in case of type II superconductivity. The 
investigation of nonuniform (vortex) states can be made with the help of gradient 
terms in the free energy JH] which take into account the spatial variations of the order 
parameter field if). This task is beyond our present consideration. Rather we investigate 
the basic problem about the possible stable uniform (Meissner) superconducting phases 
which may coexist with uniform ferromagnetic order. For this aim the free energy (1) 
is quite convenient. 

In case of a strong easy axis type of magnetic anisotropy, as is in UGe 2 [Hj, the 
overall complexity of mean-field analysis of the free energy f(ip,M) can be avoided 
by performing an "Ising-like" description: M = (0,0, Ai). Further, because of the 
equivalence of the "up" and "down" physical states (±M) the thermodynamic analysis 
can be performed within the "gauge" Ai > 0. But this stage of consideration can also 
be achieved without the help of crystal anisotropy arguments. When the magnetic 
order has a continuous symmetry one may take advantage of the symmetry of the total 
free energy f(ip,M) and avoid the consideration of equivalent thermodynamic states 
that occur as a result of the respective symmetry breaking at the phase transition point 
but have no effect on thermodynamics of the system. In the isotropic system one may 
again choose a gauge, in which the magnetization vector has the same direction as z-axis 
(\M\ = M z = Ai) and this will not influence the generality of thermodynamic analysis. 
Here we shall prefer the alternative description within which the ferromagnetic state 
may appear through two equivalent "up" and "down" domains with magnetization Ai 
and — Ai, respectively. 
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For our aims we use notations in which the number of parameters is reduced. Intro- 
ducing the parameter 

b=(b s + u s + v a ) > (2) 

we redefine the order parameters and the other parameters in the following way: 

!fj = &V4^. = 0. e e, ; M = b i/* M ; (3) 
a s a f u s v s 

r= 7r ^jrr w= t- v= t- 



7o 5 



7 = — -771 > 7i 



6 i/ 26 V4' - {bbf y/2- 

Having in mind our approximation of uniform ip and M and the notations (2) - (3), 
the free energy density f(if>, M) can be written in the form 

/(^,M) = r0 2 + i0 4 + 270 1 2 Msin(^-^ 1 )+ 7l 2 M 2 + tM 2 + iM 4 (4) 
-2w [<^sin 2 (0 2 - Ox) + <^sin 2 (0i - 3 ) + <^sin 2 (0 2 - 3 )] 

In this free energy the order parameters ip and M are defined per unit volume. 

In contrast to the situation in superconducting compounds, for the case of neutron 
stars, the crystal field anisotropy represented by the v s — terms in (1) - (4) can be ig- 
nored, and for this reason we shall discuss the case v s = 0. We assume that Tf > T s . 
This is the case when the superconductivity is triggered by the magnetic order. Besides 
we shall discuss the stable phases in the temperature region T > T s . The case Tf < T s 
may also present interest for neutron stars and, hence, it needs a special investiga- 
tion. As mentioned in Ref. ^H], the case T s ~ Tf allows for a quite simple analytical 
treatment. All these cases may be of interest to the description of ferromagnetic su- 
perconductivity in stellar objects whereas in condensed matter only cases of Tf ^> T s 
have been observed so far. 

Our consideration is performed within the framework of the standard mean-field anal- 
ysis P3]. The stable phases correspond to global minima of the GL energy (1). The 
equilibrium phase transition line separating two phases is defined by the thermody- 
namic states, where the respective GL free energies are equal. 

3. Phases and phase diagram 

The calculations show that for temperatures T > T s , i.e., for r > 0, we have three 
stable phases. Two of them are quite simple: the normal (N-) phase (ip = M = 0) 
with existence and stability domains given by t > and r > 0, and the ferromagnetic 
phase (FM) given by ip — and M 2 = —t which has the existence condition t < 0, and 
a stability domain defined by the inequalities r > jit and 



r>7it + 7v /= ^- (5) 
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Figure 1: Phase diagram in the (t, r) plane for 7 = 1.2, 71 = 0.8, and w = 0.4. 

The third stable phase is a phase of coexistence of superconductivity and ferromag- 
netism (hereafter referred to as FS). It is given the following equations: 

01 = 02 = -y| , 03 = , (6) 

2 = (± 7 M - r - 7l M 2 ) , (7) 
(1 - 7l 2 )M 3 ± ^77iM 2 +(t-\- j!r] M ± ^ = , (8) 



and 

(6 2 -e 1 ) = Tl + 27rk J (9) 

(fc = 0, ±1, ...). The upper sign in Eqs. (7) - (9) corresponds to a domain in which 
sin (#2 — 9\) = — 1 and the lower sign corresponds to a second domain which may be 
referred to as FS*; in the latter, sin(# 2 — 61) = 1. These two domains are equivalent 
and describe the same ordering. We shall focus on the upper sign in (7) - (9), i.e. on 
FS. 

The phase diagram (t, r) is outlined in Figs. 1 and 2 for different values of the anisotropy 
parameter w. The phase transition lines for w > and w < shown in Figs. 1 
and 2, respectively, have qualitatively the same shape as the phase transition lines 
corresponding to w = ^H] but there are essential quantitative differences between 
these cases. We shall discuss them in the next section. 
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Figure 2: Phase diagram in the (t, r) plane for 7 = 1.2, 71 = 0.8, and w 
various lines are explained in the text. 



-2. The 



Note, that the phase diagrams in Figs. 1 and 2 exhibit two types of phase transitions. 
The dashed curves indicate second order phase transitions of type N-FM, FM-FS and 
N-FS, whereas the curve AC and the straight line BC indicate the first order phase 
transitions N-FS and FM-FS, respectively. The points A and B are tricritical whereas 
C is a triple point, where N, FM and FS coexist [TT]. The negative values of the 
parameter r are restricted by r(T) > r(T — 0). Having in mind this condition as well 
as the shape of the phase diagram (Figs. 1 and 2) we easily conclude that both FM 
and FS ground states (at T = 0) are possible in systems described by the model (1). 
This is the case for the itinerant magnets, mentioned above, and one may speculate 
that this situation may occur in neutron stars, too. Whether FM and FS ground state 
will occur depends on the particular values of the material parameters (T s , Tf, a s , ...). 

The final aim of the phase diagram investigation is the outline of the (T, P) diagram. 
Important conclusions about the shape of the (T, P) diagram can be made from the 
form of the (t, r) diagram without an additional information about the values of the 
relevant material parameters (a s , a/, ...) and their dependence on the pressure P. For 
example, the equilibrium temperature Tps of the phase transitions to FS phase varies 
with the variation of the system parameters (a s ,af, ...) from values which are much 
higher than the characteristic temperature T s up to zero temperature. 

4. Anisotropy effects 

Our analysis demonstrates that when the anisotropy of the Cooper pairs is taken in 
consideration, there will be not drastic changes in the shape the phase diagram for 



6 



r > and the order of the respective phase transitions. Of course, there will be some 
changes in the size of the phase domains and the formulae for the thermodynamic 
quantities. This is seen from Figs. 1 and 2 which are shown for the first time in the 
present report. Besides, it is readily seen from Figs. 1 and 2 that the temperature 
domain of first order phase transitions and the temperature domain of stability of 
FS above T s essentially vary with the variations of the anisotropy parameter w. The 
parameter w will also insert changes in the values of the thermodynamic quantities 
like the magnetic susceptibility and the entropy and specific heat jumps at the phase 
transition points 21. 

Besides, and this seems to be the main anisotropy effect, the w- and v -terms in the free 
energy lead to a stabilization of the order along the main crystal directions which, in 
other words, means that the degeneration of the possible ground states is considerably 
reduced. This means also a smaller number of marginally stable states. 

The dimensionless anisotropy parameter w can be either positive or negative depending 
on the sign of u s . Obviously when u s > 0, the parameter w will be positive too 
(0 < w < 1). We shall illustrate the influence of Cooper-pair anisotropy in this case. 
The order parameters (M, <pj, 9j) are given by Eqs. (6), (9), 

, 2 ±7M - r - 71M 2 



(1-w) 

and 



>0, (10) 



3 

w - 7l 2 )M 3 ± -771M 2 + 



7 2 

t(l-w) - — -7ir 



M±^ = 0, (11) 



where the meaning of the upper and lower sign is the same as explained just below 
Eq. (9). We consider the FS domain corresponding to the upper sign in the Eq. (10) 
and (11). The stability conditions for FS read, 

(2 - w)jM - r - 71M 2 . . 

^ '- >0, 12) 

1 — w 



7M - wr - w^M 2 > 



(13) 



and 



1 — w 



3(1 



w 



7i)M 2 + 377iM + t(l 



w\ 



r 

7ir 

2 n 



> . 



(14) 



For M 7^ (7/271) we can express the function r(M) defined by Eq. (10), substitute the 
obtained expression for r(M) in the existence and stability conditions (10)-(14) and 
do the analysis in the same way as for w = [TH] - The most substantial qualitative 
difference between the cases w > and w < is that for w < the stability of FS is 
limited for r < 0. This is seen from Fig. 2 where FS is stable above the straight dotted 
line for r < and t < 0. This includes into consideration also purely superconducting 
(Meissner) phases as ground states. 



5. Final remarks 



7 



We have done an investigation of the M-trigger effect in unconventional ferromagnetic 
superconductors. This effect due to the Mipi ^-coupling term in the GL free energy 
consists of bringing into existence of superconductivity in a domain of the phase di- 
agram of the system that is entirely in the region of existence of the ferromagnetic 
phase. This form of coexistence of unconventional superconductivity and ferromag- 
netic order is possible for temperatures above and below the critical temperature T s , 
which corresponds to the standard phase transition of second order from normal to 
Meissner phase - usual uniform superconductivity in a zero external magnetic field, 
which appears outside the domain of existence of ferromagnetic order. Our investiga- 
tion has been mainly intended to clarify the thermodynamic behaviour at temperatures 
T s < T < Tf, where the superconductivity cannot appear without the mechanism of 
M-triggering. We have described the possible ordered phases (FM and FS) in this most 
interesting temperature interval. 

The Cooper pair and crystal anisotropics have also been investigated and their main 
effects on the thermodynamics of the triggered phase of coexistence have been es- 
tablished. In discussions of concrete real material one should take into account the 
respective crystal symmetry but the variation of the essential thermodynamic prop- 
erties with the change of the type of this symmetry is not substantial when the low 
symmetry and low order (in both M and ip) 7-term is present in the free energy. 

Below the superconducting critical temperature T s a variety of pure superconducting 
and mixed phases of coexistence of superconductivity and ferromagnetism exists and 
the thermodynamic behavior at these relatively low temperatures is more complex than 
in known cases of improper ferroelectrics; see. e.g., Ref. [22] • The case Tf < T s also 
needs a special investigation. 

Our results are referred to the possible uniform superconducting and ferromagnetic 
states. Vortex and other nonuniform phases need a separate study. 

More experimental information about the values of the material parameters (a a , a/, ...) 
included in the free energy (1) is required in order to outline the thermodynamic 
behavior and the phase diagram in terms of thermodynamic parameters T and P. In 
particular, a reliable knowledge about the dependence of the parameters a s and aj on 
the pressure P, the value of the characteristic temperature T s and the ratio a s /aj at 
zero temperature are of primary interest. 

The phenomeno logical GL model (1) is quite general and can be reliably used in con- 
siderations of a possible coexistance of ferromagnetism and unconventional supercon- 
ductivity in stellar objects. Recent investigations [21 E] of superconductivity in neutron 
stars can be related with the present consideration. 
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